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We review basic computational techniques for simulations of various magnetic properties of solids. 
Several applications to compute magnetic anisotropy energy, spin wave spectra, magnetic suscepti- 
bilities and temperature dependent magnetisations for a number of real systems are presented for 
illustrative purposes. 

PACS numbers: 



This review covers main techniques and their appUcations developed in the past to calculate properties of magnetic 
systems by the methods of electronic structure theory of solids. The basic tool which was used in connections with 
these developments is a spin dependent version of density functional theorji which will be reviewed in Section II in its 
most general non-collinear form and including so called LDA+U methodSi. One of the oldest application of this method 
is the problem of magnetic anisotropy energy (MAE). Due to its technological relevance as well as the smallness of 
MAE, such calculations represent a real challenge to the theory and we review these efforts in Section III. In Section 
IV we address the problem of computing non-collinear spin alignments which is solved using an elegant spin spiral 
approach^. This allows us to simulate adiabatic spin dynamics of real magnets with relatively cheap computational 
effort, and some of the calculations will be discussed and compared to experiments. Another development which has 
been undertaken in the past is the direct calculation of exchange integrals by utilizing a formula of linear response 
theoryi. This is reviewed in Section V. The spin waves, their dispersions and lifetimes as well as other spin fluctuations 
formally appear in the structure of dynamical spin susceptibility. The developments of electronic structure based on 
linear response approaches which access this quantity is reviewed in Section VI. The effect of electronic correlations 
are not small in the real materials and static mean field treatment done with local spin density functional or LSDA+U 
method may not be adequate for a whole range of solids. Bringing the effects of dynamical correlations to compute 
magnetic properties of real systems has been recently undertaken by utilizing the dynamical mean field theory^. This 
approach will be reviewed in Section VII together with some of its most recent applications. The reader is also refereed 
to several other reviews^ on a similar subject as well as a few excellent books^. 



Here we describe a basic tool for computational studies of magnetic phenomena, a spin dependent version of the 
density functional theory^. We use general non-collinear notations and also treat the effects of spin -orbit coupling. 
We utilize local spin density approximation (LSDA) to the exchange-correlation functional and also show how the 
corrections due to strong electronic correlation effects appear in the functional in its simplest Hartree-Fock form 
known as the LDA+U method-^. In section V we will also describe a more general first -principle method based on 
dynamical mean field theorji^ to include the correlation effects in a more rigorous manner. 

To describe a magnetic solid we consider a system of fermions under an external potential Vext and an external 
magnetic field 'Rext- It is useful to introduce the notion of the Kohn-Sham potential 14// (r) and the Kohn-Sham 
magnetic field Be//(r). When spin-orbit coupling is present, the intra-atomic magnetization m(r) is not coUinear, 
and the solid may choose different atom dependent quantization axes which makes magnetic moments pointing in 
different directions. Therefore, the magnetization must be treated as a general vector field, which also realizes non- 
collinear intra-atomic nature of this quantity. Such general magnetization scheme has been recently discussed^^. In 
the non-collinear version of LSDA approach the total energy functional E is considered as a functional of two variables 
the charge density p{t) and magnetization density m(r) 



I. INTRODUCTION 



II. SPIN DENSITY FUNCTIONALS IN NON-COLLINEAR FORM 
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where {||} denotes averaging over spin degrees of freedom only, the spin angular momentum operator is expressed in 
terms of Pauli matrices s = o"/2, / is 2 x 2 unit matrix, g is the gyromagnetic ratio which is for electrons equal to 2, 
and ■0kj" (i") ^I'S the non-interacting Kohn-Sham particles are formally described by the two-component spinor wave 
functions 
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which define the charge and non-coUincar magnetization densities of the electrons. 

LSDA+U method^ introduces additional variable "occupancy spin density matrix" hah 
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which represents the correlated part of electron density. To build up occupancy matrix one introduces a set of localized 
orbitals (/>Q(r), associated with correlated electrons. Then 

=E/(^k,)(<.Va)(0hlV'Lf ) (5) 

The occupancy matrix becomes non-diagonal in spin space when spin-orbit coupling is taken into account. We include 
spin-orbit coupling in a variational way as suggested by Anderseni^. 

One introduces a Lagrange multipliers matrix ISVah to enforce ((SJ. The LSDA+U total energy functional is given 
by: 



i?[p(r),m(r),n,,]=E/(ek,)ek,-y' K//(r)p(r)dr + ^ Be//(r)m(r)dr- ^ ^ ACh'^n^,-' 

ah era' 

+ j V,,t{v)p{v)dv ~ j Be,*(r)m(r)dr + i j drdr' ^^'"^^^'"'^ 



|r - r' 

[p, m] + E^'"''-' [hab] - E'^S"'' [hab] ■ (6) 
The energies e\^j are determined by Pauli-like Kohn-Sham matrix equation: 

ah 

Here 1 and s are one-electron orbital and spin angular momentum operator, respectively, ^(r) determines the strength 
of spin-orbit coupling and in practice is determinedi^ by radial derivative of the I — component of the Kohn-Sham 
potential inside an atomic sphere: 

. l^kiW^ ,8) 

dr 

m] is the LSDA exchange correlation energy. When magnetization is present, the exchange-correlation 
energy functional is assumed to be dependent on density and absolute value of the magnetization: 

El:^^%M = j dre,c[p(r),m(r)]p(r) 

+ j drfxc[p{Y),ni{v)] m(r), (9) 

where TO(r) = |m(r)|. 
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jTjModei j^^^j jg ^ contribution from the Coulomb energy in the shell of correlated electrons 
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This expression is nothing else as a Hartree-Fock average of the original expression for the Coulomb interaction 

]^'^^{ah\vc\cd)cl„cl^,Cdcr'Cccr (11) 

trcr' ahcd 

where the Coulomb interaction t!c(r — r') has to take into account the effects of screening by conduction electrons. The 
matrix {ac\vc\bd) = Uatcd is standardly expressed via the Slater integrals which for d-electrons are three constants 
^(0)^ ^(2)^ and i^^^-* considered as the external parameters of the method. Its determination can for example be done 
using atomic spectral data, constrained density functional theory calculations or taken from spectroscopic experiments. 

Since part of this Coulomb energy is already taken into account in LDA functional, the double-counting part 
denoted E^g'^'^^ [fiab] must be subtracted. Frequently used form^ for E^q'^'^'' [fiab] is 

E^'S'^' inab] = lunin - 1) - ij[nT(fiT _ 1) + n^(n^ _ 1)] (12) 



where 

Tl = 

{21- 

wirT) E((«^i^i«^) - («H (14) 



U = --r^T.(-b\]-\ab) (13) 
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and where nf = n^^ and n — -\- . 

The Kohn-Sham potential Ve//(r) and Kohn-Sham magnetic field Be//(r) are obtained by extremizing the ftmc- 
tional with respect to p{r) and m(r): 

..„M^ (15, 

Extremizing with respect to n^^ yield the correction to the potential ^V^^ 



AV;r = — -, -^^^ (17) 



ab ab 

Then the spin diagonal elements for potential correction are given by: 

= Y.^-bcdnJ-'^ +Y,iUabcd~Uadcb)Kd 

cd cd 

-SabU{n - i) + <5ah J(n"" - (18) 

and spin off-diagonal elements are given by 

^VaV^ = -T.Uadcbn-r (19) 
cd 

The off-diagonal elements of the potential correction only present when spin-orbit coupling is included, hence a 
relativistic effect. 

This completes the description of the method and we now turn out to several applications of it. 
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III. MAGNETIC ANISOTROPY OF FERROMAGNETS 

One of the oldest problems which has been addressed in the past is the calculation of the magneto-crystalline 
anisotropy energy (MAE)iii^ii^iiS4ii of magnetic materials. The MAE is defined as the difference of total energies 
with the orientations of magnetization pointing in different, e.g., (001) and (111), crystalline axis. The difference is 
not zero because of spin-orbit effect, which couples the magnetization to the lattice, and determines the direction of 
magnetization, called the easy axis. 

Being a ground state property, the MAE should be accessible in principle via spin density functional theory described 
above, and, despite the primary difficulty related to the smallness of MAE (~ 1 /iteF/atom), great efforts to compute 
the quantity with advanced total energy methods combined with the development of faster computers, have seen 
success in predicting its correct orders of magnitudes However, the correct easy axis of Ni has not been 

predicted by the LSDA and a great amount of work has been done to understand what is the difficulty. These include 
(i) scaling spin-orbit coupling in order to enlarge its effect on the MA E^^i^" , (ii) calculating torque to avoid comparing 
large numbers of energy^°, (iii) studying the effects of the second Hund's rule in the orbital polarization theory ■^i, 
(iv) analyzing possible changes in the position of the Fermi level by changing the number of valence electrons—, (v) 
using the state tracking method—, and (iv) real space approach^. 

It was recently suggested that the deficiency of the calculations is lying in improper treatment of the correlation 
effects and the intra-atomic repulsion U and exchange J should be taken into account. It is important to perform the 
calculations for fixed values of magnetic moments which themselves show some dependency on U and J as studied 
previouslj»2^. Since the pure LSDA result {U=0, and J=0) reproduces the experimental values for magnetic moments 
in both Fe and Ni fairly well, the U — J parameter space should be scanned and the path of U and J values which 
hold the theoretical moments close to the experiment was extracted. 
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FIG. 1: The magneto-crystalline anisotropy energy MAE = _E(111) — iJ(OOl) for Ni as function of U. The experimental MAE 
is marked by arrow (—2.8 fieV). The values of exchange parameter J for every value of U axe chosen to hold the magnetic 
moment of 0.61 fig 

The effect of correlations was found to be crucial in predicting the correct axis of Ni. Fig.^shows the results of this 
calculated MAE as a function of Coulomb parameter U. Walking along the path of parameters U and J which hold 
the magnetic moment to 0.6 fig the MAE first increases to 60 fieV {U = 0.5 eV, J — 0.3 eV) and then decreases. 
While decreasing it makes a rather flat area from U = 1.4 eV, J = 0.9 eV to U = 1.7 eV, J = 1.1 eV where MAE is 
positive and around 10 fieV. After the flat area, the MAE changes from the wrong easy axis to the correct easy axis. 
The correct magnetic anisotropy is predicted at [/ = 1.9 eV and J = 1.2 eV. The change from the wrong easy axis 
to the correct easy axis occurs over the range of SU ~ 0.2eV, which is of the order of spin-orbit coupling constant 
(- O.leV"). 

For Fe (see Fig. O , the MAE was calculated along the path of U and J values which fixes the magnetic moment to 
2.2 fig. At U = eV and J = eV, the MAE is 0.5 fj,eV. The correct MAE with the correct direction of magnetic 
moment is predicted atU = 1.2 eV and J = 0.8 eV. It is remarkable that the values of U and J necessary to reproduce 
the correct magnetic anisotropy energy are very close to the values which are needed to describe photoemission spectra 
of these materials2&. 

Recently, there has been a lot of experimental and theoretical studies devoted to CrOa^^'^^'^Si^Si^. This compound 
has unusual half-metallic nature, it is a metal in one spin channel and an insulator in the other. Cr02 demonstrates 
ferromagnetic ordering with magnetic moment 2fif^ per Cr atom and Curie Temperature Tc — 390K . A number of 
published worksSLSSiSi^fliSi addressed its band structure. 

The results of the MAE computations within LSDA have been studied recently. Total energy calculations for 
three different directions [001], [010] and [102]. [001] axis indicate that easy magnetization axis within LSDA is 
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FIG. 2: The magneto-crystalline anisotropy energy MAE = i5(lll) — _E(001) for Fe as function of U. The experimental MAE is 
marked by arrow (1.4 fieV). The values of exchange parameter J for every value of U are chosen to hold the magnetic moment 
of 
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FIG. 3: The magneto-crystalline anisotropy energies MAE for Cr02 as functions of U. The experimental value of MAE 
-E[010] — -^[OOl] — 15.6iieV per cell is shown by arrow. 



consistent with latest thin film experiments^^iSi^. The numerical values of MAE however exceed the experimental 
one approximately two times. 

In order to figure out the infiuence of intra-atomic repulsion U on the magnetic anisotropy, calculations have been 
also performed for finite value of U changing it from to 6eV (J — Q.87eV has been kept constant). These results are 
presented in Fig. O MAE is decreasing rapidly starting from LDA value « 68^6^ per cell and changes its sign around 
U « 0.8eV. This leads to switching of correct easy magnetization axis [001] to the wrong one, namely [102]. The 
biggest experimental value of MAE reported in the literature is 15.6/iey per cell^^. The calculated MAE approaches 
this value around point U = O.GeV. To summarize, the LSDA+U approach with U « O.GeV and J = Q.87eV 
adequately describes the magneto-crystalline anisotropy of Cr02. 

We can conclude that the calculations of magnetic anisotropy energies in materials are only in its beginning stage of 
development .It is by now clear that the effects of Coulomb correlations must be taken seriously into account but the 
smallness of MAE and its sensitivity to the values of U used in the calculations prompts us that a predictive power of 
the method is yet to be achieved. More work is clearly required to develop robust algorithms calculating this highly 
important property of ferromagnets. 



IV. SPIN SPIRAL METHOD AND FROZEN MAGNON CALCULATIONS 



While ferromagnetic spin alignments are relatively simple to explore by the methods of spin density functional 
theory, the situation gets formally much more complicated when non-coUinearity occurs. The simplest example of 
non-collinearity is the antiferromagnetism which is the ordering with wave- vector q corresponding to a zone-boundary 
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point of the Brillouin zone. This requires a choice of doubled unit cell so that in the new lattice this vector q becomes 
one of the reciprocal lattice vectors. In a more general sense we can consider an appearance of a spin spiral which 
gives orientation of different magnetic moments M.r+fl for atoms of the sublattice r + R (r is the basis vector, R is 
the primitive translation) in the form 



(20) 



where Mt is the length of the atomic moment, and 0t,0t ^^e the angles describing its orientation in the unit cell 
with R = 0. As we see, in the systems with non-collinear ordering the situation gets immediately very complicated 
as for arbitrary q the size of the unit cell becomes prohibitively large. Note that here we do not assume that the size 
of the moment is varied when going from one cell to another. The latter is another type of so called incommensurate 
magnetism, the most known example of which is Cr. 

Here we will briefly review a very elegant way to solve the problem of non-collinear magnetism where the size 
of the moment is kept constant which is known as the spin spiral method developed by Sandratski. The problem 
when the size itself varies is more difficult one and the use of linear response theory will be discussed later of this 
review. The formal point is to note that non-collinear magnetic state destroys the periodicity of the original lattice 
and therefore the original Bloch theorem applied to the Kohn-Sham states '0ki('') ^lo longer available. However, the 
Bloch theorem in the present form reflects the translational symmetry properties of the lattice only it does not take 
into account the symmetry properties of the spin subsystem. A generalized symmetry treatment of the Hamiltonian 
for the non-collinear magnet can be developed. Let us define group operations acting on a spinor state '0(r). These 
are (i) the translational operator Tr so that TjiX = r + R, (ii) the generalized rotation g which is a pure rotation 
a = {ug^Pg^'-^g) described by the 3x3 j = 1 Wigner matrix U^^,{ag, l3g,jg) (where ag,l3g,jg are the Euler angles) 
and possible shift by vector b for non-symmorphic group so that gr =dr + b, (iii) the rotation ^ = {a^,(3^, 7^) in the 

spin-i subspace described by the 2x2 j — Wigner matrix ?7^~y^(a^, /3j, 7^). 

The usefulness of these operations becomes immediately transparent if we explore symmetry properties of non- 
relativistic Hamiltonians -ffcTCT' (r) — {—'V'^ +Vef f)Saai +g^Q^effSrycr' of the type l(7|l with non-collinear spin alignments 
which assumes that the effective magnetic field Bg/ / inside each atom centered at r + R transforms similar to H20|l . 
The translations R combined with the spin rotations £,{qRT} — (^T,qR + 0riO) described by the unitary matrices 
[7^=1/2 igave the spin spiral structure invariant. We thus arrive to a generalized Bloch theorem: 



4,,(r + R) = e'^^W='/\er, qR + 0,, 0)4,(r) 



(21) 



Remarkably that this property allows us to restrict the consideration to a chemical unit cell and not to the supercell. 
Hence, instead of introducing the supercell whose size depends on the symmetry of wave vector q, the non-collinear 
spin spiral states can still be treated with the Hamiltonians whose size does not depend on the wavevector q. 
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FIG. 4: Comparison between calculated using frozen magnon method and experimental^® spin wave spectrum for Fe. 



A variety of different calculations has been done with the spin spiral approach^-^. Halilov at. al^SiSi have studied 
adiabatic spin waves computed using a frozen magnon method. This is analogous to the frozen phonon method 
where atoms are displaced from their equilibrium positions and the total energy is restored as a function of these 
positions. Here the calculation is simplified since there is no necessity to introduce supercells, and can be performed 
for any general spin wave- vector q. Fig. ^reproduces the comparison between calculated using frozen magnon LSDA 
approach and experimental^- spin wave spectra for Fe. 
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FIG. 5: Comparison between calculated using LDA+U method and experimental spin wave spectrum for NiO. 



As we see the obtained agreement in these calculations is very satisfactory despite the fact that LSDA theory does 
not take into account important effects of dynamical correlations. For systems like NiO, the correlation effects become 
more important and the frozen magnon calculations using LSDA+U method have been carried out. Fig. |S1 shows 
the comparison between calculated using LSDA+U method and experimental spin wave spectra for NiO. The overall 
agreement is good but important discrepancies remain which point out on the necessity to treat correlations effects 
among d-electrons beyond its simplest Hartree-Fock treatment as it is done within LSDA+U. 



V. CALCULATIONS OF EXCHANGE CONSTANTS 



While the density functional theory provide a rigorous description of the ground state properties of magnetic 
materials, the finite-temperature magnetism is estimated following a simple suggestion^, whereby constrained DFT 
at r = is used to extract exchange constants for a classical Heisenberg model, which in turn is solved using 
approximate methods (e.g. RPA, mean field) from a classical statistical mechanics of spin systems^i^Si^SiiS. The recent 
implementation of this approach gives reasonable value of transition temperature for iron but not for nickel*^ . The 
analysis of exchange parameters for different classes of magnetic materials such as dilute magnetic semiconductors^-, 
molecular magnets-*^, colossal magnetoresistance perovskitea^, transition metal alloys*'^, hard magnetic materials 
such as PtCo^ gives a useful information for magnetic simulations. 

In this section we outline the general strategy for estimations of exchange interactions in solids. The most reliable 
way to consider spin excitations of itinerant electron magnets in the framework of the spin density functional approach 
is the use of frequency dependent magnetic susceptibility within the linear response theory^^ii^iiS. 



6m = xS'Bext = XoSBeff, (22) 

where Xo is a "bare" DFT and x is an enhanced susceptibility. In coUinear magnetic structures there are no coupling 
between the longitudinal and transverse components and for the transverse spin susceptibility we have the following 
equation: 

X+-(r,r',c.) =Xf["(r,r',c.) + J dr"xo+"(r, r", c.)/(r")x+-(r", r', c.) (23) 

where / = Bxc/m is an exchange-correlation "Hund's rule" interaction. This random-phase-approximation (RPA) 
like equation is formally exact in the adiabatic time dependent version of density functional theory (TD-DFT)^. The 
bare susceptibility has the following form: 

X+- (r, r', ..) ^ ^ ,/_'l~^_^t €l (r)^.i (r')^pT (^0 (24) 
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where ip^^ and e^o- are eigenstates and eigenvalues for the Kohn-Sham quasiparticles. We can rewrite the equation 
for the transverse susceptibihty as 

X+- ^{m + A^{uj- /,,a) (25) 

where 

A(r, r', = ^ -Al-J^4,* (r)^^^ (r) V [i',,ir')Vr.^ (r') - (r') V^^^ (r')] (26) 

Spin wave excitations can be separated from the Stoner continuum (e.g., paramagnons) only in the adiabatic 
approximation, which means the replacement A(r,r',a;) by A(r, r',0) in Ea. H25|l . Otherwise one should just find the 
poles of the total susceptibility, and the whole concept of "exchange interactions" is not uniquely defined. Nevertheless, 
formally we can introduce the effective exchange interactions via the quantities 

r!(r,r',c^) =4,A(r,r',w). (27) 

In the static limit one can show that 

n{v, r', 0) = ^ J(r, r', 0) - B,.,{r)5{r ~ r') (28) 
m[r) 

where an expression for frequency dependent exchange interactions has the following form 

J(r, r', ..) = ^ (r)i?.c(r) V-.i (rW., (r')i?.c(r')^^T (^0 (29) 

The later coincides with the exchange integrals^ if we neglect the w - dependence. Since Bxc ~ rn we have J ^ m? 
and the expression l(?7jl vanishes in non-magnetic case, as it should be. Note that the static susceptibility x^" [uj — 0) 
can be represented in the following form 

X+- (0) = m {n-^ - B-^) ~ mfl (30) 

which is equivalent to the result of Rcf.^^ 

n = ?i(i-B-}nY^ (31) 



for the renormalized exchange interaction if one define them in terms of inverse static susceptibility^^ . Note, that the 
magnon frequencies are just eigenstates of the operator Vl (0) which exactly corresponds to the expression from the 
unrenormalized exchange interactione^ . Note that for the long- wavelength limit q — + this result turns out to be 
exact which proves the above statement about the stiffness constant D: in the framework of the local approximation 
it is rigorous. 

Let us compare the Eg 13 II with the Heisenberg-like dispersion law in the spherical approximation for the spin- wave 
spectrum in terms of static transverse susceptibility*^ : 

= W - Xo"'] ■ (32) 

Now we will use an expression J = x~^which gives the well-known connection between the inverse static susceptibility 
and the parameter $2 in the Landau theory of phase transitions. To analyze the relationship between the dispersion 
law and that commonly used in the DFT, oj^ = ml [xo — Xq] we assume that the ratio 

Aq = (Xq - Xo) Xo"' ~ ^'Jrnl (33) 

is small in the long- wave approximation. Then, by expanding Eg 13 21 over the parameter Aq, we obtain the desired 
result: 

Wq = TO (xq^ - Xo^) = mXo^Aq (1 - Aq)~^ 

- mxo' (xo - Xq) Xo ' = (4 - -/q) (34) 
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where the susceptibihty has a matrix structure and 



Jo {I 



ixJ 



(36) 



is the matrix of the exchange parameter in the local (long- wave) approximation. We take into account, that in the 
ferromagnetic state Xo^ = I. A form of Ea. H36() is the widely accepted definition of the exchange coupling parameter. 
The usual assumption of weak enhancement J„ — Jg 



Jq — Jq is valid due to the well-known result Xq ' 



J„ 



As a consequence, the spectrum of elementary excitations is not affected by exchange-correlation enhancement effects 
in linear response regime. So, the definition H3t)|) is directly related to several very strong approximations: rigid spin, 
smallness of spin wave dispersion compared to the effective exchange splitting (atomic limit). 

Below we will briefly mention calculations of the exchange parameters and adiabatic spin wave spectra using the 
multiple scattering technique in the local density approximation with a linear mufRn-tin orbital technique (LMTO) in 
the atomic sphere approximationii. Three ferromagnetic systems with entirely different degree of localization of the 
local moment were considered: Gd, Ni and Fe^. An important observation is that in Gd (highly localized moments, 
small spin wave dispersion) the approximation based on the assumption Aq << f is completely fulfilled (A at the 
zone boundary is less then O.Of), whereas in Ni (moderately localized moments, large spin wave dispersion) it is 
not valid at all (A ~ 0.6 at q —Y), so that the corresponding matrix estimates using Ea. (|34|l demonstrate a strong 
enhancement of spin wave spectra at larger q. Fe is an intermediate case where the maximum of A is 0.30. 

Practical expression for the exchange coupling has been derived^ using multiple scattering method: 
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47r 
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00 
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(37) 



where T" is a path operator. This expression properly takes into account both exchange effects related to the Stoner 
splitting and the different energy dispersion for different spins in band magnets. 

Analysis of the effective exchange coupling between atoms after Eq. H37|l indicated that in Fe and Ni the main 
contribution came from the renormalization of the first nearest neighbor exchange, so that in bcc Fe it is increased 
from Jgj^=16.6 meV to Joi=f9.4 meV, whereas in fee Ni Jq^—2.7 meV to Joi=8.3 meV. Such a clear difference in 
the results indicates that the removal of the long-wave approximation can serve as an indicator of the degree of 
localization (parameter (xq — Xo) Xo^ above) in different metallic magnets. In Fe and Ni, for instance, it at least 
partially explains why long wavelength mean field estimates predict such a small Tc in FM Ni: 300-350K in Refi^, 
340 K in present calculations, while experimental result is 630 K. First of all, long-wave approximation is suitable 
for such a 'localized' system as Fe, and the corresponding change in the nearest neighbor Jqi is relatively small 
(correspondingly the increase of Tc is expected to be small). Ferromagnetic Ni represents a rather itinerant system 
and any local approach (long wave approximation in particular) might produce a large error. Also in Ni at Tc the 
energy associated with 'creation' of moment is comparable with the energy of its rotation and longitudinal fluctuations 
should be taken into account. 

A large increase in Jqi for Ni may indicate that traditional mean-field approaches (based on the absence of short- 
range order) are not applicable for the itinerant systems in general, predicting very high Tc (above lOOOK). Such 
a number is not consistent with the local density approach because the latter does not allow to have Tc in fee Ni 
larger then 500-540K (600-640K using gradient corrections) and this number can be considered only as an indicator 
of non-applicability of mean field approach. 



VI. DYNAMICAL SPIN SUSCEPTIBILITY CALCULATIONS 



Local spin density functional theory and the LDA+U method have been applied with success to study full 
wave-vector and frequency dependent spin susceptibility x of solids. Its knowledge accessible directly via neutron- 
scattering measurements is important due to significant influence of spin fluctuations to many physical properties 
and phenomena*, such, e.g., as the electronic specific heat, electrical and thermal resistivity, suppression of supercon- 
ductivity for singlet spin pairing, etc . As we already discussed, in magnetically ordered materials, transverse spin 
fluctuations are spin waves whose energies and lifetimes are seen in the structure of transverse susceptibility. 

Large efforts in the past were put on the development of methods for ab initio calculations of the dynamical spin 
susceptibility based either on the random-phase-approximation (RPA) decoupling of the Bethe-Salpeter equation^—, 
or within density functional formalism.^"* . Quantitative estimates of x with realistic energy bands, wave functions, and 
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self-consistently screened electron-electron matrix elements appeared in the literature more than 3 decades ago^SiSi^. 
More recently, a computationally efficient time-dependent generalization of the Sternheimer approach^ originally 
developed for ab initio calculations of phonon dispersions, electron-phonon interactions and transport properties of 
transition-metal materials has been proposed^. 

This method is advantageous since it allows us to avoid problems connected to the summations over higher energy 
states and inversions of large matrices. It considers the response of electrons to a small external magnetic field 

SBextirt) = ,5be^('i+^)''e"^*e-''l*l + c.c (38) 

applied to a solid. Here Sh — dh^^e^ shows a polarization of the field (/i runs over x, y, z or over —1, 0, 1), wave 
vector q lies in the first Brillouin zone, G is a reciprocal lattice vector, and 77 is an infinitesimal positive quantity. 
If the unperturbed system is described by charge density p(r) and by magnetization m(r), the main problem is to 
find self-consistently first-order changes 6p{rt) and Sm{rt) = (5m^(rt)e'' induced by the field SHextiJ^t). If the 
polarization Sh in H38|l is fixed to a particular /zth direction, and Sm.{rt) is calculated afterwards, a /ith column of the 
spin susceptibility matrix Xi^^(r, q + G, will be found. 

Employment of time-dependent (TD) version of density functional theory (DFT)-^*^- to find the quantities Sp(rt) 
and (5m(ri) can be useful because in order to find the dynamical response within TD DFT, only the knowledge of 
these unperturbed Kohn-Sham states (both occupied and unoccupied) is required; no knowledge of real excitation 
spectra (both energies and lifetimes) is necessary. This is the main advantage of such approach. Unfortunately, 
within TD DFT, an accurate approximation to the kernel /a;c(r, r',a;) describing dynamical exchange-correlation 
effects is unknown while some progress is currently been madeS. In the past, the so called adiabatic local density 
approximation (ALDA)S and a generalized gradient approximation^ (GGA) are adopted to treat /a;c(r, r', w). The 
addition of the local U corrections has been addressed in the worliS. 

One can develop a variational linear-response formulation. It has already been known that static charge and spin 
susceptibilities appeared as second-order changes in the total energy due to applied external fields can be calculated 
in a variational way. This was demonstrated long time agc^^ on the example of magnetic response, and, recently^^'^^ , 
in the problem of lattice dynamics which is an example of charge response. The proof is directly related to a powerful 
"2n+ 1" theorem of perturbation theory and stationarity property for the total energy itself^^. Any (2n+ l)th change 
in the total energy Etot involves finding only (n)th order changes in one-electron wave functions i/'i, and corresponding 
changes in the charge density as well as in the magnetization. Any (2n)th change in Etot is then variational with 
respect to the (n)th-order changes in tp^. 

A time-dependent generalization of these results was addressed^, where the action 5 as a functional of p{rt) and 
m(rt) is considered within TD DFTiSiSi. These functions are expressed via Kohn-Sham spinor orbitals i/j-^^j {rt) 
satisfying TD Schrodinger's equation. Therefore, S as the stationary functional of t/jj.^ (rt) is considered in practice. 
When the external field is small, the perturbed wave function is represented as ip^j (r) e^''^^^^ + Sip^{rt) and the 
first-order changes Sip^^j (rt) define the induced charge density as well as the magnetization: 

Sp = + (39) 

ki 



Sm = .gAis^ ({'5V^kj |§IV^kj} + Wkj l§l'54j}) (40) 



In order to find (5'0kj" (''^)i a time-dependent analog of the "2n+l" theorem is now introduced. Any (2n+l)th change 
in the action functional S involves finding only (n)th order changes in the TD functions V'kj (i"^)i and corresponding 
changes in charge density as well as in the magnetization. Any (2rt)th change in S is then variational with respect to 
the (n)th-order changes in tp^^j (rt). The proof is the same as for the static case®'^ if the stationarity property of S and 
the standard TD perturbation theory are exploited. For important case n = 2, this theorem makes the second-order 
change S*^^^ in the action variational with respect to the first-order changes S'tpi (rt) . If the perturbation has the form 
(|38|l . S'^^^ is directly related to the real diagonal part of the dynamical spin susceptibility i?e[x,y^(q+G', q+G, cj)]g'=g 
, thus allowing its variational estimate^. 

The problem is now reduced to find S^'^^ as a functional of S'tpi^j (rt) and to minimize it. This will bring an equation 

for S-tpY^j (rt). Any change in the action functional can be established by straightforward varying S of TD DFTi^SiSi 
with respect to the perturbation H38|l . This is analogous to what is done in the static DFT to derive, for example, 
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the dynamical matrix^. 5^^^ is found to be 



y SpSVeff - J Sin{5B,ff + 5B,,t) (41) 

where the unperturbed 2x2 Hamihonian matrix H = (— V^ + Vg//)/— 5/i^sBe//+<^ls. SV^ff and <^Be// are their first- 
order changes induced by the perturbation 1)38(1 which involve the Hartree (for SVeff) and the exchange-correlation 
contributions expressed via (5/9 and i5m in the standard manner^. 

The differential equation for Sip^^j (rt) derived from the stationarity condition of H41|l is given by 

{H - id^ji)di^^j + {sv.ffi - ^^B(TSB,ff)lp^,^ = o (42) 



Applications to transverse spin fluctuations in Fe and Ni as well as calculations of paramagnetic response in Cr 
and Pd have demonstrated an efficiency of the approach. Experimental evidence of an " optical" spin-wave branch 
for Ni^^ and its absence for Fe2& was correctly described . The dynamical susceptibility was calculated ab initio for 
paramagnetic Cr, a highly interesting material due to its incommensurate antiferromagnetism-^. The calculation 
predicted a wave vector of the spin density wave (SDW), and clarified the role of Fermi-surface nesting. Strong 
long-wavelength spin fluctuations of Pd were evident from these and earlier— theoretical studies. 




We reproduce some of the calculations at Fig. Eland Fig. [7| which show calculated Im[x{(l,uj)] for ferromag- 
netic Ni and paramagnetic bcc Cr. In particular, for Cr, a remarkable structure is clearly seen for the q's near 
{0,0,xsDW ~0.9)27r/a, where the susceptibility is mostly enhanced at low frequencies. This predicts Cr to be an 
incommensurate antiferromagnet (experimentally, xsdw=0-95). 

A scheme for making ab initio calculations of the dynamic paramagnetic spin susceptibilities of solids at finite tem- 
peratures was recently described^ where incommensurate and commensurate antiferromagnetic spin fluctuations in 
paramagnetic Cr and compositionally disordered CrgsVs and CrgsRes alloys were studied together with the connection 
with the nesting of their Fermi surfaces. 

To conclude, the developed approach is able to describe known spin-fluctuational spectra of some real materials 
which demonstrated its efficiency for practical ab initio calculations, however, more elaborate approximations to the 
dynamical exchange and correlation are clearly required in order to account for the observed discrepancies. The next 
section discussed a more general framework to deal with the correlation effects. 
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FIG. 7: Calculated 7m[x+_ (q, t^)] (Ry^^) for Cr. 



VII. DYNAMICAL MEAN FIELD APPROACH AND THE REGIME OF STRONG CORRELATIONS 



Magnetic materials range from very weak having a small magnetization to strong ones which exhibit a saturated 
magnetization close to the atomic value. As we have seen in the previous sections, frequently, weak or itinerant 
magnets are well described by spin density wave theory, where spin fluctuations are localized in a small region of 
momentum space. Quantitatively they are well described by LSDA. The transition from ordered to paramagnetic 
state is driven by amplitude fluctuations. In strong magnets, there is a separation of time scales, h/t is the time scale 
for an electron to hop from site to site with hopping integral t, which is much shorter than h/ J, the time scale for the 
moment to flip in the paramagnetic state. The spin fluctuations are localized in real space and the transition to the 
paramagnetic state is driven by orientation fluctuations of the spin. The exchange splitting J is much larger than the 
critical temperature. 

Obtaining a quantitative theory of magnetic materials valid both in weak and in strong coupling, above and below 
the critical temperature has been a theoretical challenge for many years. It has been particularly difficult to describe 
the regime above Tc in strong magnets when the moments are well formed but their orientation fluctuates. This 
problem, for example, arises in magnetic insulators above their ordering temperature, when this ordering temperature 
is small compared to electronic scales, a situation that arises in transition metal monoxides (NiO and MnO) and led 
to the concept of a Mott insulator. In these materials the insulating gap is much larger than the Neel temperature. 
Above the ordering temperature, we have a collection of atoms with an open shell interacting via superexchange. This 
is a local moment regime which cannot be accessed easily with traditional electronic structure methods. 

Two important approaches were designed to access the disordered local moment (DLM) regime. One 
approach^ starts form a Hubbard like Hamiltonian and introduces spin fluctuations via the Hubbard-Stratonovich 
transformationSiffliiiiZ^ which is then evaluated using a static coherent potential approximation (CPA) and improve- 
ments of this technique. A dynamical CPA^'^ was developed by Kakehashi^^. A second approach begins with solutions 
of the Kohn-Sham equations of a constrained LDA approximation in which the local moments point in random di- 
rections, and averages over their orientation using the KKR-CPA approach2Si2&. The average of the Kohn-Sham 
Green functions then can be taken as the first approximation to the true Green functions, and information about 
angle resolved photoemission spectra can be extracted^LS. There are approaches that are based on a picture where 
there is no short range order to large degree. The opposite point of view, where the spin fluctuations far away form 
the critical temperature are still relatively long ranged was put forward in the fluctuation local band picturrfMSsii. 

Description of the behavior near the critical point requires renormalization group methods, and the low temperature 
treatment of this problem is still a subject of intensive research^. There is also a large literature on describing magnetic 
metals using more standard many-body mcthodaSM^iSSiSSiSi. 

Dynamical mean field theory (DMFT) can be used to improve treatment to include dynamical fluctuations beyond 
static approximations like LSDA or LSDA+U. This is needed to get access to disordered local moment regime and 
beyond. Notice that single site DMFT includes some degree of short range correlations. Cluster methods can be used 
to go beyond the single site DMFT to improve the description of short range order on the quasiparticle spectrum. 
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DMFT also allows us to incorporate the effects of the electron-electron interaction on the electronic degrees of freedom. 
This is relatively important in metallic systems such as Fe and Ni and absolutely essential to obtain the Mott-Hubbard 
gap in transition metal monoxides. 

Realistic implementation of DMFT within electronic structure framework^*'^^ can be viewed as a spectral density 
functional theory^Si^iiSSiS^. The central quantity of this method is the local Green function G;oc(r, r', w), which in the 
non-coUinear version should be considered as the matrix in spin variables 

'^loc ^ ( I (43) 

The local Green function is only a part of the full one-electron Green function restricted which is restricted within a 
certain cluster area. Due to translational invariance of the Green function on the original lattice given by primitive 
translations {R}, i.e. G(r + R, r' + R, w) = G(r,r',a;), it is always sufficient to consider r lying within a primitive 
unit cell J^c positioned at R = 0. Thus, r' travels within some area fi/oc centered at R = 0. We therefore define the 
local Green function to be the exact Green function G(r,r',z) within a given cluster lljoc and zero outside. In other 
words, 

Gioc{r, r', u) = G(r, r', u)eioc{r, r') (44) 

where the theta function is a unity when vector r G fie, r' £ ^loc and zero otherwise. A functional theory where the 
total free energy of the system Tsdf considers the local Green function as a variable was developed recentljiSfliSiiSSiS^. 
It is more powerful than the ordinary DFT since it accesses both the energetics and local excitational spectra of real 
materials. The range of locality is in principle our choice and the theory has a correct scaling in predicting full k- 
and uj dependent spectrum of excitations when this range is expanded till infinity. If the range is collapsed to a single 
site, only the density of states is predicted. 

A numerically tractable approach was developed completely similar to the density functional theory with the 
introduction of energy-dependent analog of Kohn-Sham orbitals. They are very helpful to write down the kinetic 
energy portion of the functional similar to as DFT considers the density functional as the functional of Kohn-Sham 
wave functions. It is first useful to introduce the notion of an local self-energy operator 

= [Ve^t{v)+VH{v)]6{v~v')6,„,+MZ\vy,uj) 

which has the same range of locality as the local Green function. Second, an auxiliary Green function Q'^'^ (r, r', iui) 
connected to our new "interacting Kohn-Sham" particles so that it is defined in the entire space by the relationship 

[g-r^'(r,r',c.) = [Go-r'^'(r,r',^) - A^^„-;(r,r',a;). 

where Gq"^ (r, r', iuj) is the non-interacting Green function, which is given by 

[Go w) = [(w + /i + - Kx*(r))(5,,- - ^1 • ^aa' - g^iB'Bexts^.'W " r') (45) 

Since Q is a, functional of Gioc, it is very useful to view the spectral density functional T sdf as a functional of Q: 

Tsdf[Q] =T:T\ng -Ty[G^^ ^ g-^]g + '^sdf[Gioc] (46) 

where the first two contribution represent the kinetic contribution and the contribution from the external potential 
and magnetic field. The unknown interaction part of the free energy ^sdf[Gioc] is the functional of Gioc- If the 
Hartree term is explicitly extracted, this functional can be represented as 

^SDF [Gioc] = Eh [p] + ^foF [Gioc] (47) 



To facilitate the calculation, the auxiliary Green function Q'^'^ (r,r',iw) can be represented in terms of generalized 
energy-dependent one-particle states 

(r, r', lu) = y "^k." ^ ^"^k,^ V ; ^^g^ 
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which obey the one-particle Dyson equation for the "right" eigenvector i^^j^\'r) 

^{[-V^ + K.t(r) + VH{r)]S,,, + g^isBeZ/s,,, + ^1 ■ s..' Wkf'^W + (49) 

a' 

^ J dv'MZ'{v,v'Mi^^:\v') = E^,^i,^l\v) (50) 

cr' 

and similar equation for the "left" eigenvector V'kj^'' ('^) when it appears on the left. 

Evaluation of the free energy requires writing down the precise functional form for '^sdf[Gioc\- Unfortunately, it is 
a problem because the free energy F = E — TS, where E is the total energy and 5 is the entropy, and the evaluation 
of the entropy requires the evaluation of the energy as a function of temperature and an additional integration over 
it. However, the total energy formula can be written by utilizing the Galitzky-Migdal expression for the interaction 
energy^. It is given by: 

E = e'-o^ 5kj-i?k,c. -TJ2T.I dvdv'Mlff{T, r', iu)g''' {v' , r, iu:) + 

+ J drVeMr)pir)- J drBe^t{r)mir)+EH[p] + lTj2Yl j drdv' {v,v' ,iuj)g^' ^ [v\v,ilo) (51) 



where 

ffkjo; = -. — ; — ^—^ (52) 

We have cast the notation of spectral density theory in a form similar to DFT. The function gkjcj is the Green function 
in the orthogonal left/right representation which plays a role of a "frequency dependent occupation number". 
Both the charge density and non-coUinear magnetization can be found from the integral over the k-space: 

- J^EEE.9'^-<t^W<t^(r)e-°" (53) 



m(r) =5Msr^^.gk,.^s..,V^^]:)(r)^g:''(r)e-°" (54) 

The central procedure which is yet to be determined is the construction of the exchange-correlation part of the 
self-energy. It is formally a variational derivative of the exchange correlation part of the free energy functional 

M^c'iryM = = ^^^f^M,v') (55) 

The spectral density functional theory, where an exact functional of certain local quantities is constructed in the 
spirit of Ref. |93 uses effective self-energies which are local by construction. This property can be exploited to 
find good approximations to the interaction energy functional. For example, if it is a priori known that the real 
electronic self-energy is local in a certain portion of the Hilbert space, a good approximation is the corresponding 
local dynamical mean field theory obtained for example by a restriction or truncation of the full Baym-Kadanoff 
functional to local quantities in the spirit of Ref. I^The local DMFT approximates the functional ^sdf[Gioc] by 
the sum of all two-particle diagrams evaluated with G/oc and the bare Coulomb interaction vc ■ In other words, the 
functional dependence of the interaction part $ bk [G] in the Baym-Kadanoff functional for which the diagrammatic 
rules exist is now restricted by G/oc and is used as an approximation to ^sdf[Gioc], i-e. ^sdf[Gioc] = ^bk[Gioc]- 
Obviously that the variational derivative of such restricted functional will generate the self-energy confined in the 
same area as the local Green function itself. 

Remarkably the summation over all local diagrams can be performed exactly via introduction of an auxiliary 
quantum impurity model subjected to a self-consistency condition^. If this impurity is considered as a cluster G, the 
cellular DMFT (C-DMFT) can be used which breaks the translational invariance of the lattice to obtain accurate 
estimates of the self energies. 
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Various methods such as LDA+U2i, LDA+DMFT~ and local GW22i2^ which appeared recently for realistic calcu- 
lations of properties of strongly correlated materials can be naturally understood within spectral density functional 
theory. Let us, for example, explore the idea of expressing the energy as the density functional. Local density ap- 
proximation prompts us that a large portion of the exchange-correlation part <&a;c[p] can be found easily. Indeed, 
the charge density is known to be accurately obtained by the LDA. Why not think of LDA as the most primitive 
impurity solver, which generates manifestly local self-energy with localization radius collapsed to a single r point? It 
is tempting to represent ^sdf[Gioc] = Eh[p] + E^^^[p] + ^Gioc] - ^dc[Gioc], where the new functional ^sdf[Gioc] 
needs in fact to take care of those electrons which are strongly correlated and heavy, thus badly described by LDA. 
Conceptually, that means that the solution of the cluster impurity model for the light electrons is approximated by 
LDA and does not need a frequency resolution for their self-energies. 

Such introduced LDA-I-DMFT approximation considers both the density and the local Green function as the 
parameters of the spectral density functional . A further approximation is made to accelerate the solution of a 
single-site impurity model: the functional dependence comes from the subblock of the correlated electrons only. If 
localized orbital representation {Xa} is utilized, a subspace of the heavy electrons {Xa} can be identified. Thus, the 
approximation can be written as ^SDF[Gioc,ab{iuj)] , where Gioc,ab{i^) is the heavy block of the local Green function. 

Unfortunately, the LDA has no diagrammatic representation, and it is difficult to separate the contributions from 
the light and heavy electrons. The Ef^l^^[p\ is a non-linear functional and it already includes the contribution to 
the energy from all orbitals in some average form. Therefore we need to take care of a non-trivial double counting, 
encoded in the functional ^dc[Gioc\- The precise form of the double counting is related to the approximation imposed 
for ^[Gioc] and can be borrowed from the LDA+U method^. Note that both LDA+U and LDA+DMFT methods 
require separations of the electrons onto light and heavy which makes them basis-set dependent. 

Iron and nickel were studied in Refs. [96 97. The values U = 2.3 (3.0) eV for Fe (Ni) and the same value of the inter- 
atomic exchange, J = 0.9 eV for both Fe and Ni were used, as came out from the constrained LDA calculations^SiSSiSS. 
These parameters are consistent with those of many earlier studies and resulted in a good description of the physical 
properties of Fe and Ni. In Ref. 97 the general form of the double counting correction V^'-^ — \Trc^J^a{Q) was taken. 
Notice that because of the different self energies in the eg and t2g blocks the DMFT Fermi surface does not coincide 
with the LDA Fermi surface. 

The impurity model was solved by QMC in Ref. |93 and by the FLEX scheme in Ref. llOOl It is clear that nickel 
is more itinerant than iron (the spin-spin autocorrelation decays faster), which has longer lived spin fluctuations. On 
the other hand, the one-particle density of states of iron resembles very much the LSDA density of states while the 
DOS of nickel below Tc, has additional features which are not present in the LSDA spectrai2i4i224i2^: the presence 
of a famous 6 eV satellite, the 30% narrowing of the occupied part of d-band and the 50% decrease of exchange 
spittoons compared to the LDA results. Note that the satellite in Ni has substantially more spin-up contributions in 
agreement with photoemission spectraifi^. The exchange splitting of the d-band depends very weakly on temperature 
from T = 0.6Tc to T = 0.9rc- Correlation effects in Fe are less pronounced than in Ni, due to its large spin splitting 
and the characteristic bcc structural dip in the density of states for the spin-down states near the Fermi level, which 
reduces the density of states for particle hole excitations. 

The uniform spin susceptibility in the paramagnetic state, Xq=Q — dM / dH , was extracted from the QMC simulations 
by measuring the induced magnetic moment in a small external magnetic field. It includes the polarization of the 
impurity Weiss field by the external field^. The dynamical mean field results account for the Curie- Weiss law which 
is observed experimentally in Fe and Ni. As the temperature increases above Tc, the atomic character of the system 
is partially restored resulting in an atomic like susceptibility with an effective moment: 

Xg=o = 2.{T%cy ^^^^ 

The temperature dependence of the ordered magnetic moment below the Curie temperature and the inverse of the 
uniform susceptibility above the Curie point are plotted in Fig. |S| together with the corresponding experimental 
data for iron and nickel^i. The LDA-I-DMFT calculation describes the magnetization curve and the slope of the 
high-temperature Curie- Weiss susceptibility remarkably well. The calculated values of high-temperature magnetic 
moments extracted from the uniform spin susceptibility are Me^y = 3.09 (1.50)//^ for Fe (Ni), in good agreement with 
the experimental data Me// — 3.13 (1.62)/!^ for Fe (Ni)i2^. 

The Curie temperatures of Fe and Ni were estimated from the disappearance of spin polarization in the self- 
consistent solution of the DMFT problem and from the Curie- Weiss law in Ea. (|56|l . The estimates for Tc — 1900 
(700) K are in reasonable agreement with experimental values of 1043 (631) K for Fe (Ni) respectively^^i., considering 
the single-site nature of the DMFT approach, which is not able to capture the reduction of Tc due to long-wavelength 
spin waves. These effects are governed by the spin-wave stiffness. Since the ratio of the spin-wave stiffness D to Tc, 
Tc/a?D, is nearly a factor of 3 larger for Fe than for NiiSi (a is the lattice constant), we expect the Tc in DMFT to 
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FIG. 8; Temperature dependence of ordered moment and the inverse ferromagnetic susceptibility for Fe (open square) and Ni 
(open circle) compared with experimental results for Fe (square) and Ni (circle) (from Ref. ? ). The calculated moments were 
normalized to the LDA ground state magnetization (2.2 for Fe and 0.6 /ig for Ni). 

be much higher than the observed Curie temperature in Fe than in Ni. Quantitative calculations demonstrating the 
sizeable reduction of Tc due to spin waves in Fe in the framework of a Heisenberg model were performed in Ref. l4ll 
This physics whereby the long wavelength fluctuations renormalize the critical temperature would be reintroduced in 
the DMFT by the use E-DMFT. Alternatively, the reduction of the critical temperature due to spatial fluctuations 
can be investigated with cluster DMFT methods. 

Within the dynamical mean field theory one can also compute the local spin susceptibility defined by 



is easily computed within the DMFT-QMC. Its behavior as a function of temperature gives a very intuitive picture 
of the degree of correlations in the system. In a weakly correlated regime we expect Ea. (|57f) to be nearly temperature 
independent, while in a strongly correlated regime we expect a leading Curie- Weiss behavior at high temperatures 
Xiocai — M?oc/("^-^ "I" const.) where ^i^^ is an effective local magnetic moment. In the Heisenberg model with spin S, 
l^foc ~ '^'('5'+ 1)5^ and for the well-defined local magnetic moments (e.g., for rare-earth magnets) this quantity should 
be temperature independent. For the itinerant electron magnets, fii^^ is temperature dependent due to a variety of 
competing many-body effects such as Kondo screening, the induction of local magnetic moment by temperature* and 
thermal fluctuations which disorder the momentsiS^. All these effects are included in the DMFT calculations. 

The comparison of the values of the local and the q — Q susceptibility gives a crude measure of the degree of 
short-range order which is present above Tc- As expected, the moments extracted from the local susceptibility 
Ea. (|57|l are a bit smaller (2.8 fig for iron and 1.3 fig for nickel) than those extracted from the uniform magnetic 
susceptibility. This reflects the small degree of the short-range correlations which remain well above To^SL. The 
high-temperature LDA-I-DMFT clearly shows the presence of a local moment above Tc- This moment is correlated 
with the presence of high-energy features (of the order of the Coulomb energies) in the photoemission. This is 
also true below Tc, where the spin dependence of the spectra is more pronounced for the satellite region in Nickel 
than for that of the quasiparticle bands near the Fermi level. This can explain the apparent discrepancies between 
different experimental determinations of the high-temperature magnetic splittingsiSSiiSSiiSSiiiS as being the results 
of probing different energy regions. The resonant photoemission experiments^^^ reflect the presence of local-moment 
polarization in the high-energy spectrum above the Curie temperature in Nickel, while the low-energy angle resolved 
photoemission investigationgWO, results in non-magnetic bands near the Fermi level. This is exactly the DMFT view 
on the electronic structure of transition metals above Tq. Fluctuating moments and atomic-like conflgurations are 
large at short times, which results in correlation effects in the high-energy spectra such as spin-multiplet splittings. 




(57) 







where S = i a a' Cma^aa'Cma' IS singlc-sitc Spin Operator. It differs from the q — susceptibility by the absence 
of spin polarization in the Weiss field of the impurity model. Ea. H57|) cannot be probed directly in experiments but it 
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The moment is reduced at longer time scales, corresponding to a more band-like, less correlated electronic structure 
near the Fermi level. 

NiO and MnO represent two classical Mott-Hubbard systems. Both materials are insulators with the energy gap 
of a few eV regardless whether they are antiferro- or paramagnetic. The spin dependent LSDA theory strongly 
underestimates the energy gap in the ordered phase. This can be corrected by the use of the LDA+U method. Both 
theories however fail completely to describe the local moment regime reflecting a general drawback of band theory to 
reproduce atomic limit. Therefore the real challenge is to describe the paramagnetic insulating state where the self- 
energy effects are crucial both for the electronic structure and for recovering the correct phonon dispersions in these 
materials. The DMFT calculations have been performediii by taking into account correlations among d-electrons. 
In the regime of large U adequate for both for NiO and MnO in the paramagnetic phase the correlations were treated 
within the well-known Hubbard I approximation. 

The calculated densities of states using the LDA+DMFT method for the paramagnetic state of NiO and MnOiii 
have revealed the presence of both lower and upper Hubbard subbands. These were found in agreement with the 
LDA+U calculations of Anisimo\«^ which have been performed for the ordered states of these oxides. Clearly, spin 
integrated spectral functions do not show an appreciable dependence with temperature and look similar below and 
above phase transition point. 

The progress in the electronic structure calculation with DMFT allowed to apply the method to more complicated 
systems such as plutoniumiiSiii^. The temperature dependence of atomic volume in Pu is anomalousii^iii^. It shows 
an enormous volume expansion between a and S phases which is about 25%. Within the S phase, the metal shows 
negative thermal expansion. Transition between S and higher-temperature e phase occurs with a 5% volume collapse. 
Also, Pu shows anomalous resistivity behavior^^^ characteristic for the heavy fermion systems, but neither of its 
phases is magnetic. The susceptibility is small and relatively temperature independent. The photoemission^^^ shows 
a strong narrow Kondo-like peak at the Fermi level consistent with large values of the linear specific heat coefficient. 

Density functional based LDA and GGA calculations describe the properties of Pu incorrectly. They predict 
magnetic orderingii^. They underestimateiiS equilibrium volume of the 6 and e phase by as much as 30%, the 
largest discrepancy known in LDA, which usually predicts the volume of solids within a few % accuracy even for such 
correlated systems as high temperature superconductors. 

To address these questions several approaches have been developed. The LDA+U method was applied to S Pui2S. 
It is able to produce the correct volume of the S phase for values of the parameter J7~4 eV consistent with atomic 
spectral data and constrained density functional calculations. Similar calculation has been performed by a so-called 
orbitally ordered density functional methodiSi. However, both methods predict Pu to be magnetic, which is not 
observed experimentally. The LDA+U method is unable to predict the correct excitation spectrum. Also, to recover 
the a phase within LDA+U the parameter U has to be set to zero which is inconsistent with its transport properties 
and with microscopic calculations of this parameter. Another approach proposediS^ in the past is the constrained LDA 
approach in which some of the 5f electrons, are treated as core, while the remaining are allowed to participate in band 
formation. Results of the self-interaction-corrected LDA calculations have been reportedi^^, as well as qualitative 
discussion of the bonding nature across the actinides series has been giveoi^i. 

Dynamical mean-field calculations have been recently applied with success to this material. They were performed 
in paramagnetic state and do not impose magnetic ordering. DMFT was able to predict the appearance of a strong 
quasiparticle peak near the Fermi level which exists in the both phases. Also, the lower and upper Hubbard bands 
were clearly distinguished. The width of the quasiparticle peak in the a phase is found to be larger by 30 per cent 
compared to the width in the S phase. This indicates that the low-temperature phase is more metallic, i.e. it has larger 
spectral weight in the quasiparticle peak and smaller weight in the Hubbard bands. Recent advances have allowed 
the experimental determination of these spectra, and these calculations were consistent with these measurementsiii. 

VIII. CONCLUSION 

In conclusion, electronic structure methods have been applied with success to study magnetic properties of solids. 
Many good results have been obtained for weakly correlated itinerant systems. New algorithms for inclusion of many- 
body effects are beginning to attack magnetic phenomena in strongly correlated materials. These successes have 
generated a new synergy between the electronic structure and many-body communities, which are jointly developing 
new and powerful electronic structure tools for virtual material exploration. 
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